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A simple explicit example of a Roberts-type dynamo is given in which the a-efFect of mean-field 
electrodynamics exists in spite of point -wise vanishing kinetic helicity of the fluid flow. In this way 
it is shown that a-effect dynamos do not necessarily require non-zero kinetic helicity. A mean-field 
theory of Roberts-type dynamos is established within the framework of the second-order correlation 
approximation. In addition numerical solutions of the original dynamo equations are given, that 
are independent of any approximation of that kind. Both theory and numerical results demonstrate 
the possibility of dynamo action in the absence of kinetic helicity. 
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I. INTRODUCTION 

The essential breakthrough in the understanding of the 
origin of the large-scale magnetic fields of cosmic objects 
came with the development of mean-field dynamo theory. 
A central component of this theory is the a-effect, that 
is, a mean electromotive force with a component parallel 
or antiparallel to the mean magnetic field in a turbulcntly 
moving electrically conducting fluid. The a-effect, which 
occurs naturally in inhomogeneous turbulence on a rotat- 
ing body, is a crucial element of the dynamo mechanisms 
proposed and widely accepted for the Sun and planets, 
for other stellar objects and even for galaxies; see, e.g.. 

In the majority of investigations this effect has been 
merely calculated in the so-called second-order correla- 
tion approximation (SOCA), or first-order smoothing ap- 
proximation (FOSA), which ignores all contributions of 
higher than second order in the turbulent part of the fluid 
velocity. Moreover in many cases attention has been fo- 
cussed on the high-conductivity limit only, which can 
be roughly characterized by short correlation times of 
the turbulent motion in comparison to the magnetic- 
field decay time for a turbulent eddy. Under these cir- 
cumstances the a-effect is closely connected with some 
average over the kinetic helicity of the turbulent mo- 
tion, that is, of M • a;, where u means the turbulent 
part of the fluid velocity, u) the corresponding vortic- 
ity, u; = V X u; see, e.g., P, SSflHl- More precisely, 
the coefficient a for isotropic turbulence turns out to be 
equal to — i u{x, t) ■ u){x, t — t) At, often expressed 
in the form — ■ uj Tc with equal arguments of u and u) 
and some appropriate time Tc. Here and in what follows 
overbars indicate averages. For anisotropic turbulence 
the trace of the a-tensor is just equal to three times this 
value of a. 

These findings have been sometimes overinterpreted in 
the sense that mean-field dynamos, or even dynamos at 
all, might not work without kinetic helicity of the fluid 



motion, that is, if it • a; vanishes. There exist however 
several counter-examples. 

Firstly, the mean-field dynamo theory offers also dy- 
namo mechanisms without a-effect, for which the av- 
erage of the kinetic helicity may well be equal to zero. 
For instance, a combination of the so-called fi x J- 
effect, which may occur even in the case of homoge- 
neous turbulence in a rotating body, with shear asso- 
ciated with differential rotation, can act as a dynamo 
IS S i, S E, [H [H, [H, [II- The possibility of a dy- 
namo due to turbulence influenced by large-scale shear 
and the shear flow itself [1^ [3 [13] is still under debate 
j8, 19, 20, 21]. 

Secondly, in dynamo theory beyond the mean-field 
concept a number of examples of dynamos without ki- 
netic helicity are known. The dynamo proposed by 
Herzenberg [22| , usually considered as the first existence 
proof for homogeneous dynamos at all, works without ki- 
netic helicity. Likewise in the dynamo models of Gailitis 
working with an axisymmetric meridional circulation in 
cylindrical or spherical geometry [H, [2^ there is no 
kinetic helicity. We mention here further the dynamo in 
a layer with hexagonal convection cells proposed by Zhe- 
ligovsky and Galloway [26|], in which the kinetic helicity 
is zero everywhere. 

Thirdly, it is known since the studies by Kazantsev [2^] 
that small-scale dynamos may work without kinetic he- 
licity. This fact has meanwhile been confirmed by many 
other investigations; see, e.g., [1]. 

Let us return to the mean-field concept but restrict 
ourselves, for the sake of simplicity, to homogeneous and 
statistically steady turbulence. We stay first with the 
second-order correlation approximation but relax the re- 
striction to the high-conductivity limit. As is known 
from early studies, e.g. [l|, [H, [23], the crucial parame- 
ter for the a-effect is then no longer the mean kinetic 
helicity u ■ uj. The a-effect is in general determined by 
the function /i($,t) = u{x,t) ■ uj{x + + t) or, what 
is equivalent, by the kinetic helicity spectrum, defined 



as its Fourier transform with respect to ^ and r. In the 
high-conductivity hmit this brings us back to the above- 
mentioned results. In the low conductivity limit, that is, 
large correlation times in comparison to the magnetic- 
field decay time for a turbulent eddy, the a-effect is 
closely connected with the quantity u ■ xp, where xp is the 
vector potential of u, that is 'V x if) — u with ^ -ip — 0. 
For isotropic turbulence the coefficient a is then equal to 
— ■^u-tp, where rj means the magnetic diffusivity of the 
fluid [l|, 0, [H, [2^ . In the anisotropic case the trace 
of the a-tensor is again equal to three times this value 
of a. In general u ■ uj does not vanish in this limit but 
it is without interest for the a-effect. Both u ■ uj and 
u ■ xp are quantities, which, if non-zero, indicate devia- 
tions of the turbulence from reflectional symmetry. In 
the second-order correlation approximation the fl x J- 
effect, too, is completely determined by the kinetic helic- 
ity spectrum [T^ . 

Beyond the second-order correlation approximation 
the a-effect is no longer determined by the kinetic he- 
licity spectrum alone. In an example given by Gilbert 
et al. [30| an a-effect (with zero trace of the a-tensor) 
occurs in a non-steady flow even with zero kinetic helic- 
ity spectrum. Note that small-scale dynamos may also 
work with zero kinetic helicity spectrum, but they do not 
produce large scale fields. 

Interesting simple dynamo models have been proposed 
by G. O. Roberts (sil . [3^ . In his second paper [s^l he 
considered dynamos due to steady fluid flows which are 
periodic in two Cartesian coordinates, say x and y, but 
independent of the third one, z. When speaking of a 
"Roberts dynamo" in what follows we refer always to the 
first flow pattern envisaged there [equation (5.1), figure 
1]. This dynamo played a central role in designing the 
Karlsruhe dynamo experiment [H, HI, [H, [s^ . It can be 
easily interpreted within the mean-field concept. In that 
sense it occurs as an a-effect dynamo with an anisotropic 
a-effect. The flow pattern proposed by Roberts, and in a 
sense realized in the Karlsruhe experiment, shows non- 
zero kinetic helicity. In this paper we want to demon- 
strate that this kind of dynamo works also with a slightly 
modified flow pattern in which the kinetic helicity is ex- 
actly equal to zero everywhere (but the kinetic helicity 
spectrum remains non-zero). 

In section II we present a simple mean-field theory of 
Roberts dynamos at the level of the second-order correla- 
tion approximation. In section III we report on numerical 
results that apply independently of this approximation. 
Finally, in section IV some conclusions are discussed. 

II. A MODIFIED ROBERTS DYNAMO 

A. Starting point 

We focus our attention now on the Roberts dynamo in 
the above sense and refer again to a Cartesian coordinate 
system {x,y,z). The fluid is considered as incompress- 
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FIG. 1: Tlie Roberts flow pattern 

ible. Therefore its velocity, u, has to satisfy V • m = 0, 
and it is assumed to be a sum of two parts, one propor- 
tional to e X Vx and the other to xe, where e is the unit 
vector in z direction, and x — sin^Trx/a) sin{TTy /a). The 
corresponding flow pattern is depicted in figure [1] We 
denote the sections defined by na < x < {n + l)a and 
ma < X < (m + l)a with integer n and m as "cells". In 
that sense the velocity u changes sign when we proceed 
from one cell to an adjacent one. A fluid flow of that kind 
acts indeed as a dynamo. Magnetic field modes with the 
same periodicity in x and y varying however with z like 
sin(fcz) or cos(fc2) may grow for sufficiently small values 
of |fc|. The most easily excitable modes for a given k pos- 
sess a part which is independent of x and y, but varies 
with z. 



B. Mean— field theory with a generalized fluid flow 

Let us proceed to more general flow patterns which 
show some of the crucial properties of the specific flow 
pattern considered so far. We assume again that the flow 
is incompressible, steady, depends periodically on x and y 
but is independent of z, and that it shows a cell structure 
as in the example above. More precisely we require that 
u changes sign if the pattern is shifted by a length a along 
the X or y axes or is rotated by 90° about the z axis. 

Let us sketch the mean-field theory of dynamos work- 
ing with fluid flows of this type (partially using ideas 
described in our paper [371] , in the following referred to 
as RB03). 

It is assumed that the magnetic field B is governed by 
the induction equation 

rjV^B + V X {ux B) - dtB = 0, V ■ B = 0, (1) 
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with the magnetic difFusivity 77 being constant. 

We define mean fields by averaging the original fields 
over a square that corresponds to four cells in the xy 
plane. More precisely we define the mean field which 
belongs to an original field -F, in any given point x = x^y 
by averaging F over the square defined hy Xq — a < x < 
xo + a and yo — a < y < yo + a. If applied to quantities 
which are periodic in x and y with a period length 2a 
this average satisfies the Reynolds rules. Clearly we have 
u = 0. 

When taking the average of equation ([T]) and denoting 
B — B hy b we obtain the mean-field induction equation 



V xe^dtB^O, V ■B = Q, 



with the mean electromotive force 



£ = u X b. 



(2) 



(3) 



For the determination of £ for a given u the equation for 
b is of interest. This equation follows from ([T]) and ([2]), 



r/V^b + V X (m X b)' - dtb = -V x {u x B) 



V b 



0, 



(4) 



where (m x 6)' stands for u x b — u x b. 

For the sake of simplicity we consider here only the 
steady case. We introduce a Fourier transformation with 
respect to z, that is, represent any function F{x,y,z) 
by F{x,y,z) — J F{x,y,k) exp{ikz) dk. If F(x,y,z) is 
real, we have F*(x,y, k) = F{x,y, — fc), where the aster- 
isk means complex conjugation. We have now 



£ = u X b . (5) 

Transforming ^ accordingly we obtain 

r){V^ - k^)b + (V + ifce) x (tt x b)' 

= -(V + ifce) X (it X B) , 
V-b + ikb, =0. (6) 

We further suppose that B is independent of x and y. 
This implies that £ is independent of x and y, too. It is 

then obvious that the connection between £ and B must 
have the form 



Si (lij B j , 



(7) 



where a^- is, like £ and B, independent of x and y, but 
depends on k. Because £*{k) = £i{—k) and B^{k) — 

Bi{~k) it has also to satisfy a*j{k) — aij{~k). 

Clearly is determined by u. Since £ does not 
change under inverting the sign of u, and a 90° rota- 
tion of the u field about the z axis changes nothing else 
than its sign, the components of the tensor dy have to 
be invariant under 90° rotations of the coordinate system 
about the z axis. We may therefore conclude that 

Oiij ^ ai{\k\) 5ij + a2{\k\) e^Cj + i asd/c]) fc ey^efe (8) 



with real ai, 02 and as. Together with ([7]) this leads to 
£z — (oi + 0'2)Bz- On the other hand, £z is equal to 
the average of u^by — Uyhx, and we may conclude from 

^ that bx and by are independent of Bz- This implies 
ai + 02 = 0. We may then write 

dy = -a^{k) {Sij - eiCj) + i /3(fc)fceijfcefe , (9) 

with two real quantities a±_ and /?, which are even func- 
tions of k. From ([7]) and (O we obtain 



£ ^ -a±{k){B- {e- B)e) -i(3kex B . 



(10) 



This relation is discussed in some detail in RB03, section 
V B. 

We now restrict our attention to the limit of small fc, 
that is, small variations of B with z. Then a± and /? 
lose their dependence on k. Denoting them by a± and P 
we conclude from (flOl) that 



£ = -c 



{B-{e-B)e) ~(3exdB/dz. (11) 



Equation ^ together with pT|) allow solutions of the 
formB = ^Boexp{ikz+Xt)] with A = ±a^k-{r]+f])k^ . 
Growing solutions are possible if 



\k\{r, + p) 



> 1 



(12) 



We note that, in agreement with the anti-dynamo theo- 
rem by Zeldovich, A vanishes if fc — » 0, that is, if B loses 
its dependence on z. 



C. Calculation of a± and f3 

For the determination of the coefficients a± and (3 we 
restrict ourselves to an approximation which corresponds 
to the second-order correlation approximation. It is de- 
fined by the neglect of (m x b)' in ([l]), or (it x b)' in ([5]). 
In view of a± and (3, which correspond to the limit of 

small k, we expand b and B in powers of k but neglect 
all contributions with higher than first powers of k. That 

.(0) ^(1) — —iO) ^(1) 

is, b = b +kb and B = B +kB , where of course 



b'"^ = b^") and b""' = From m follows then 



-(0) 



{0) 



,V^b'°' = -(^'^V)n 



9-(l) 

f]V^b ' 



—iO) ^(1) 

-ie x (m X B ) - (-B • V) m . (13) 



We put now u — 'V x xp, where if) denotes a vector po- 
tential satisfying V • = 0, further ~ 'V x ip with 
V • ■0 — 0, which implies 



(14) 
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From (flS]) with u expressed in this way we conclude 
,-(0) 



1 ,^(0) 

-(B ■V)ip 

V 



^(1) ^ 1 



. ^(0) ^(1) 

iex {'4) X B )~{B ■ . (15) 



-(0) -(1) 

Putting in the above sense also S = S + kS we find 



u X {ip X B 



(16) 



+u X {B 



V)t/,} 



We compare this with what follows from (fTU]) in this 
expansion, 



p.(0) 

t = —a I 



■^(0) ^(0) 

B {e -B ) e 



-(0) 



S = -i/3e X B 
In this way we obtain first 



B 



-(1) 



. —(1), 
(e-B ) e 



.(17) 



1 



/3 



= — {ux {g- V)xp) ■ g 

r] 

1 



7] 



u ■ ijj — [h ■ u) {h ■ tp) 



(18) 



^(0) 



where g and are unit vectors in the directions of B 

^(0) ^(0) 

[e-B ) e and ex B , respectively. Of course, a± and 
(3 cannot really depend on these directions, and we may 
average the two expressions for a± with g — (1, 0, 0) and 
g = (0, 1,0) and likewise those for f3 with h — (1,0,0) 
and h = (0, 1,0). This yields 



^M- (V X ^), 



P = -—{u-iP + {e-u)ie-xP)). (19) 



The result for a\ can also be written in the form 



1 



1 



aA_ ^ — u ■ xp ^ — ■ {V X ^p) . 

Zl] Zl] 



(20) 



In agreement with our remarks in the introduction the 
a-effect is not determined by the average of the kinetic 
helicity u ■ u) but by the related but different quantity 
u ■ ip, that is, some kind of mean helicity of the vector 
potential ^p. 




FIG. 2; A single cell of the modified Roberts flow 



the context of the Karlsruhe dynamo experiment [38| . 
To define the flow first in a single cell like those depicted 
in figure [2] we introduce a cylindrical coordinate system 
(r, if, z) with r = in the middle of this cell and z as in 
the Cartesian system used above. We put then 

Ur — everywhere 

= , Uz ^ u in r < ''I 

— , Uz = in ri < r < r2 (21) 

— LOT , Uz — euja/2 in r2 < r < 

with uj and e being constants. The full flow pattern is de- 
fined by continuation of that in the considered cell to the 
other cells such that the flow has always different signs 
in two adjacent cells. Thinking of the "central channels" 
and the "helical channels" of the Karlsruhe device we la- 
bel the flow regions r < ri and r2 < r < rj, hy C and 7i, 
respectively. Clearly finite kinetic helicity exists only in 
7i, and only for non-zero e. Note that, in contrast to the 
situation depicted in figure [1] here the kinetic helicity is 
non-negative as long as e is positive. 

For the flow defined in this way quantities like u ■ u 
ov u ■ xp depend on u and oj. In the following we express 
u by the flow rate Vc of fluid through the cross section 
< r < T'l at a given z, and w by the flow rate Vu the 
fluid through a meridional surface defined by r2 < r < 
and < z < a and a given ip. As can be easily verified 
we have then 



(22) 



with some positive quantity c depending on r2, and a. 
For u ■ ■0 it is important that the part of the vector po- 
tential xp resulting from the flow in the region C of a given 
cell docs not vanish in ?i, and likewise the part resulting 
from the flow in H, does not vanish in C. Moreover, the 
parts of xp resulting from flows in other cells have to be 
taken into account. Considering these aspects we may 
conclude that 



D. A specific flow 



U-Xp = (ciVc + £C2Vh) Vh 



(23) 



We elaborate now on the above result for a± for flow 
patterns similar to those which have been considered in 



with non-zero quantities ci and C2 depending on ri, r2, 
ra and a. 
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TABLE I: Examples of marginal dynamo states with k — 0.9 



case 


2ri/a 


2r2/a 


2r3/a 


e Vc 


Vn 


(i) 


0.5 


0.5 


1.0 


0.228 2 


0.736 


(ii) 


0.5 


0.5 


1.0 


2 


0.805 


(iii) 


0.5 


0.7 


1.0 


2 


0.965 



2.5 








2.0 






/c=1.8 ; 


1.5 








1.0 




l\ V,^^^^^^^ /c=o.9 : 


0.5 






0.0 







1 2 3 4 5 
Vc 



was chosen, for the cases (ii) and (iii) zero kinetic helicity. 
The marginal modes turned out to be non-oscillatory. 
Note that the parameters for case (i), 2ri/a = 2r2/a = 
0.5, 2r3/a = 1, e = 0.228 and k = 0.9, are realistic values 
in view of the Karlsruhe dynamo device. In the cases (i) 
and (ii) the regions C and Ti. have a common surface. In 
order to make sure that there is no effect of an artificial 
kinetic helicity due to the restricted numerical resolution 
at this surface, in case (iii) a clear separation of these two 
regions was considered. 

As long as the kinetic helicity in H does not vanish, 
that is e 7^ 0, the dynamo works even with Vc — 0, that 
is zero flow in C. In the absence of kinetic helicity, e = 0, 
the dynamos works, too, as long as Vc and therefore u ■ xjj 
are different from zero. However, the dynamo disappears 
if, in addition to e, also Vc and so u ■ ip vanish. This 
is plausible since then choices of the coordinate origin 
X = are possible such that u{x) and —u{—x) coincide. 
In this case, usually referred to as "parity-invariant" , any 
a-effect can be excluded. 

The dissipation of the mean magnetic field grows with 
K. By this reason the dynamo threshold grows with k, 
too. 



FIG. 3: Lines of marginal stability in the Vc Vh plane for the 
cases (i) (dashed lines) and (iii) (solid lines) defined in Table U 
and three values of k. 



Clearly, u ■ ip and, according to ((20|) . a± remain in 
general different from zero if the flow loses its kinetic 
helicity, that is, as u ■ w, or e, vanish. This implies the 
possibility of a-effect dynamos without kinetic helicity. 



III. NUMERICAL EXAMPLES OF q-EFFECT 
DYNAMOS WITH AND WITHOUT KINETIC 
HELICITY 

In our above calculations of a± and f3 an approxima- 
tion in the spirit of the second-order correlation approx- 
imation has been used. The results are reliable only if 
a suitably defined magnetic Reynolds number is much 
smaller than unity. This can also be expressed by re- 
quiring that the normalized flow rates Vc = Vc/jya and 
Vh = Vn/ria show this property. In addition we have 
assumed weak variations of the mean magnetic field in 
the z-direction, more precisely k = k/a 1. 

In order to confirm the existence of a-effect dynamos 
without kinetic helicity in an independent way and to 
check whether it occurs only under the mentioned condi- 
tions or over a wider range of parameters, equation ([1]) 
with the flow deflned by (PT|) has been solved numerically. 
The same numerical method as in RB03 has been used. 

In table |T] and flgure [3] some numerically determined 
dimensionless flow rates Vc and Vh are given for marginal 
dynamo states. For case (i) a non-zero kinetic helicity 



IV. CONCLUSIONS 

We have examined a modified version of the Roberts 
dynamo with a steady fiuid fiow as it was considered in 
the context of the Karlsruhe dynamo experiment. By 
contrast to what might be suggested by simple explana- 
tions of this dynamo the necessary deviation of the fiow 
from refiectional symmetry is not adequately described 
by the kinetic helicity u ■ uj oi the fluid flow or its aver- 
age u ■ ijj. Dynamo action with steady fluid flow requires 
in truth a deviation which is indicated by a non-zero 
value of the quantity u ■ ip, where ip is the vector po- 
tential of u. Even if the kinetic helicity is equal to zero 
everywhere, dynamo action occurs if only the quantity 
u ■ xp, the helicity of the vector potential of the flow, is 
unequal to zero. This fact is not surprising in light of the 
general findings of mean-field electrodynamics. 

We stress that our result does not mean that the quan- 
tity M • i/j is in general more fundamental for the a-effect 
and its dynamo action than u ■ u). In cases with steady 
fluid flow indeed u ■ xp is crucial. However, with flows 
varying rapidly in time, the crucial quantity is u ■ oj. 

Both the mean kinetic helicity u ■ uj and the quantity 
u ■ xp are determined by the helicity spectrum of the flow. 
That is, the Roberts dynamo does not work with zero 
kinetic helicity spectrum. In this sense the situation is 
different from that in the interesting example of an a- 
effect dynamo given by Gilbert et al. [soj, which works 
with a very special non-steady flow for which the kinetic 
helicity spectrum is indeed zero. 

Looking back at the Karlsruhe experiment we may 
state that a modifled (technically perhaps more com- 
plex) version without kinetic helicity of the fluid flow 
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should show the dynamo effect, too. As table |T] exem- 
plifies, however, the kinetic helicity can well reduce the 
dynamo threshold. In that sense we may modify the title 
of the paper of Gilbert et al. and say: Helicity is unnec- 
essary for Roberts dynamos - but it helps! 
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